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Abstract 

We review the formalism of primordial black holes (PBHs) produc- 
tion and show that the mass variance at horizon crossing has been 
systematically overestimated in previous studies. We derive the cor- 
rect expression. The difference is maximal at the earliest formation 
times and still very significant for PBH masses ~ 10 15 g, an accurate 
estimate requiring numerical calculations. In particular, this would 
lead to weaker constraints on the spectral index n. We then derive 
constraints on inflationary models from the fact that primordial black 
holes must not overdose the Universe. This is done both for the scale- 
free case of the power spectrum studied earlier and for the case where 
a step in the mass variance is superimposed. In the former case we 
find various constraints on n, depending on the parameters. In the 
latter case these limits can be much more strengthened, so that one 
could find from an observational limit on n a constraint on the allowed 
height of the step. 

PACS Numbers: 04.62. +v, 98.80.Cq 



1 Introduction 



Cosmology has now entered an exciting stage where observations of ever increas- 
ing accuracy can probe the remote past of our Universe. The present paradigm 
makes use of an inflationary stage of expansion in the very early Universe (see 
e.g. H, U). This solves some of the big problems of the old big bang cosmology 
in an elegant way, in particular the problems of horizon and flatness. But most 
importantly, and there lies its predictive power, it gives a possible solution to 
the crucial problem of where the primordial fluctuations leading to the observed 
large-scale structure (LSS) come from. In fact, they have their origin in the ubiq- 
uitous vacuum fluctuations. The seed of the LSS has been observed in the form of 
tiny fluctuations imprinted on the cosmological microwave background (CMB) at 
the time of decoupling. Each inflationary model makes precise predictions about 
the spectrum of its primordial fluctuations and this is how these models can be 
constrained by observations. 

It was realized already some time ago that primordial fluctuations can lead 
to the formation of primordial black holes (PBHs) (see e.g. |4j] for a review). For 
this, one needs a spectrum of primordial fluctuations whose statistics and ampli- 
tude are known. Inflationary models produce a spectrum of fluctuations also on 
scales relevant for the formation of PBHs. The statistics for PBH formation is 
usually taken as Gaussian for simplicity, which fits most of the inflationary mod- 
els. One therefore possesses an additional tool to constrain inflationary models. 
The a priori interesting thing about PBHs is that they probe scales which are 
many orders of magnitude smaller than scales probed by large-scale structure 
(LSS) surveys and CMB angular anisotropy observations. In this sense it is anal- 
ogous, even if less spectacular, to the primordial gravitational wave background 
of inflationary origin. 

It was already shown that PBH formation can put constraints on the spectral 
index of the primordial density perturbations, n < (1.23 — 1.25), see e.g. |5|, 
|6], [?|, [8]]. Basically, these constraints come either from the evaporation of black 
holes (leading to a somewhat stronger limit) or from the fact that PBHs must not 
overdose the present Universe, i.e., from the requirement Qpbh,o < 1- On mass 
scales M < 10 35 g, these constraints on n are even stronger than those coming 
from COBE. There are also constraints from the need not to overproduce PBHs 



during the preheating phase following inflation g, |10|. Such constraints assume 
of course that the primordial spectrum is scale-free, an assumption that we shall 
relax here. Hence there are two attractive aspects in the study of PBHs using 
inflationary perturbations: it expands the range of scales that can be probed, and 
it can test possible features in the spectrum on these scales not probed by CMB 
and LSS observations. 

A characteristic scale can be generated, for example, in two-fields inflation 



such as double inflation (see e.g. ITT], |12|) or in single-field inflation with broken 



scale invariance (BSI) (see e.g. (|T3|, |Ij|). Such a model was used recently in order 
to drastically cut the power on small scales in an attempt to explain the dearth of 
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dwarf galaxies [JL5|J . We want here to extend this study to the formation of PBHs. 
We are, therefore, interested in features on much smaller scales than those which 
we have considered previously, an exciting possibility which we seriously consider 
as there is no reason why the spectrum should be scale-free on all those scales 
relevant for PBH formation (though the actual occurrence of the scale-free case 
cannot be excluded). 

Our paper is organized as follows. In Sect. II we start by reviewing the basic 
formalism for PBH formation. We then point out in detail how expressions for 
the mass variance used earlier have to be corrected in order to obtain accurate 
numbers. We then consider several concrete models, beginning with the usual 
scale-free case, then extending our study to the presence of a step in the mass 
variance. In Sect. Ill we then obtain constraints on the spectral index from the 
requirement that PBHs must not overdose the Universe. For the scale-free case 
we get various constraints on n, depending on the parameters of interest, while 
for the case of a scale-free spectrum with a superimposed step the constraint on n 
depends on the height of the step. Therefore, one can obtain constraints on this 
height from the existing limits on n. Finally, Sect. IV contains our conclusions. 



2 PBH formation 
2.1 Basic formalism 

A spectrum of primordial density fluctuations can lead to the production of pri- 
mordial black holes. Let us assume for simplicity that a PBH is formed when 
the density contrast averaged over a volume of the (linear) size of the Hubble 
radius satisfies 5 min < 5 < S max , and the PBH mass, M PBH , is of the order of the 
"horizon mass" M#, the mass contained inside the Hubble volume. Usually one 
takes 5 m in = §, bmax = 1 coming from semianalytic considerations ||. Recent 
numerical calculations indicate, however, that PBHs can be formed over a much 



wider range of masses at a given formation time and that 8 m in ~ 0.7 |Tjj 17]. This 
is why we leave the limits open at this stage. For the purpose of studying the 
changes introduced by a characteristic scale in the primordial spectrum, it will 
be sufficient in the following to make the simplifying assumption that PBHs are 
formed with the mass M#. We shall, however, also comment below on the changes 
introduced by the assumption that the PBH mass is only a certain fraction of 

Each physical scale R(t) is defined by some wavenumber k and evolves with 
time according to R(t) = a(t)/k. For a given physical scale R(t), the "horizon" 
crossing time - here we do not mean "horizon" crossing during inflation, but 
after inflation - is the time when this scale reenters the Hubble radius, which 
will inevitably happen after inflation for scales that are greater than the Hubble 
radius at the end of inflation. It is at this time tk where the above condition 
for the density contrast holds. It will lead to the formation of a PBH with 
mass M PBH , being approximately equal to M H {tk). Clearly, there is a one-to-one 
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correspondence between R(tk), Mn{tk), and k. Of course we can also take this 
correspondence at any other initial time U and relate the physical quantities at 
both times ti and t k . 

Generally, if the primordial fluctuations obey Gaussian statistics, the prob- 
ability density Pr(S), where 5 is the density contrast averaged over a sphere of 
radius R, is given by 

Pn(8) = fn — ~ , p , ■ (!) 

V27T a(R) 



Here, the dispersion (mass variance) <J 2 (R) = (^(jjf^ ^ is computed using a 
top-hat window function ||], 

1 f°° 

AR) = / dA; fc , ( 2 ) 

where P(fc) is the power spectrum. It is defined by P(k) = (|<5fc| 2 ) in the case of 
finite volume. For inflation, one has to take instead the infinite-volume expression 

(^)EF(^(k-k'), (3) 

where we assume isotropy of the ensemble. From a fundamental point of view, 
the averages (...) refer to quantum expectation values; however, an effective 
quantum-to-classical transition is achieved during inflation fl8 , l9| . In the case 
of PBHs produced by inflationary perturbations, this quantum-to-classical tran- 
sition is guaranteed for the masses of interest to us (see pOf and the remarks at 
the end of section 2.2). 

The expression WthI^R) stands for the Fourier transform of the top-hat 
window function divided by the probed volume Vw = |vri? 3 , 

W TH (kR) = jt ~B^(sin kR - kR cos kR) . (4) 

[Kit I 

Therefore the probability /3(Mjj) that a region of comoving size R = k~ x has 
an averaged density contrast at horizon crossing t k in the range 8 min < 5 < S max , 
is given by 



P(M H ) = / e 2CT W dS « 1} k> e , (5) 

where cr 2 H {tk) : = <J 2 (R)\ tk , and the last approximation is valid for <5 mi „ ^> anitk), 
and (S max - 5 min ) > a H {t k ). 

Important conclusions can be drawn from (gj). Let us consider first the value of 
(3(Mh) today. We then have cr 2 H {to) ~ 1CT 8 . So clearly the probability of forming 
a black hole today is extraordinarily small. This probability can be larger in the 
primordial universe if the power is increased when we go backwards in time, but 
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usually the number will remain very small, (3{Mh) <C 1, at all times. This is due 
to the amplitude of S^^/aj^itk), on which (3(Mh) depends sensitively, see fl5|). 

The expression for (3{Mh) is usually interpreted as giving the probability 
that a PBH will be formed with a mass Mp^u > Mu{tk), i-6. greater or equal 
to the mass contained inside the Hubble volume when this scale reenters the 
Hubble radius. Strictly speaking, though, this is not true, since (^j) does not take 
into account those regions that are underdense on a scale Mh, but nevertheless 
overdense on some larger scale. In the Press-Schechter formalism (see e.g. [Q, |^]) 
this seems to be taken care of in some models by multiplying fl5|) with a factor 2. 
Fortunately, as emphasized in 0], in most cases [3{Mh) is a very rapidly falling 
function of mass, so this effect can be neglected. In this case, (3(Mh) does give 
the probability for PBH formation and thus also (at time tk) the mass fraction 
of regions that will evolve into PBHs of mass greater or equal to Mjj- This will 
also be the case for our models. 



2.2 Improved calculation of the mass variance 

As is clear from the previous subsection, one is interested in computing c 2 (R) at 
horizon crossing using the power spectrum of the primordial fluctuations P(k,t) 
of interest. Let us consider the simplest case which is usually considered when 
the power spectrum is scale-free and of the form P(k,t) = A(t) k n . Then, the 
following equation is usually used to relate <7ff(t*.) to the present value cr H (to), 
see e.g. 0, 



a 



' M H (to) ' 


n-l 

3 


~M H {t eq ~ 


M H (t eq )_ 




M H (t k )_ 



(6) 



The expression (|]) rests on the assumption that anitk) is well approximated by 
the integral (0), without window function but with a cut-off at the horizon scale. 
Consequently, the presence of the window function would be equivalent to the 
introduction of an effective cut-off at the horizon scale, leading to (recalling that 
A(t k ) oc k~ 4 ) 

a 2 H (t k ) oc k^ 1 . (7) 

However, for a scale-free primordial spectrum with n > 1, one recognizes that (0) 
is ultraviolet divergent! Hence it is clear that the integral is dominated by the 
contribution from small scales. Actually, even for 0.8 < n < 1, the small scales 
will still contribute significantly to the integral and even dominate it for n ~ 1. 
On the other hand, the quantities k 3 $ 2 (A;,t), or Sj I (k,t), defined as 



{aH) 4 k 3 
k 4 2^ 



P(k,t) 



9tt 5 



k 3 ^ 2 {k,t) 



(8) 



are time-independent on "superhorizon" scales (scales bigger than the Hubble 
radius, k < aH) in the cases that we consider here, being equal in very good 
approximation to their value at the horizon crossing time On these scales, 
these quantities behave as in (0). In addition, the constant of proportionality is 
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different in the matter and radiation dominated stages. With t kr and t km denoting 
times in the radiation and matter dominated stage, respectively, we get (cf. e.g. 
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kf^(K,t kr ) = (y) 2 ^ $2 (^^)(^ n_1 

S 2 H (kr,t kr ) = (y) 2 ^.O^)" 1 ■ (9) 

However, inspection of (||) shows that it is the quantity auitk) which is needed. 

Actually (and fortunately), there is a natural upper cut-off in &;-space for the 
power spectrum, namely k e , corresponding to the Hubble radius at the end of 
inflation t e . The lower limit can be taken zero if we assume that the number of 
e-folds during inflation amply solves the cosmological particle horizon problem. 
We thus have to consider only wavelengths bigger than k~ l . The relation between 
o\ff(^fc) and fca <3>(/c,t fc ), or dnitk), depends on the time tk and is model- dependent. 
We have in general (with k = (aH)\ tk ) 

a 2 H (t k ) = a 2 (k) 6 2 H {k,t h ) , (10) 

where 

~ he. 

a 2 (k) = / k x n+2 T 2 (kx,t k ) W^ H (x) dx . (11) 
Jo 

The transfer function T(k, t) is defined through 

P(k, t) = P(k, U) T 2 (k, t) , T(k - 0, t) - 1 . (12) 

Here, U is some initial time when all scales are outside the Hubble radius, k aH, 
and we can take ti = t e . 

The problem in evaluating a(k) comes from the evolution of the perturbations 
for scales k! inside the Hubble radius, k = (aH)\ tk < k! < k e , or equivalently 
1 < x < %■. This small scale evolution is encoded in the transfer function T(k, t) 
which is defined through relation flT2"|). For scales outside the horizon, the transfer 
function equals one. We emphasize that in flll] ), the transfer function must be 
taken at the time t k of interest, not today. Clearly, an accurate value for ( |TT1) can 
be obtained only numerically and with an explicit knowledge of T(k,t), which is 
of course model dependent j2|. 

Let us now compare our results with what is usually done using (|). The 
quantity cr^(t ) is certainly finite, even for n > 1, because one must use in (§), 
according to ([11]), the actual power spectrum whose deformation on small scales 
is encoded in the present-day transfer function T(k, to), 

poo 

a 2 H (t ) = S 2 H (t ) / x n+2 W 2 H (x) T 2 (k x,t ) dx , (13) 
Jo 
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where k = (aH) to , and Sn(t ) = <5ff(/co, t ). The transfer function today, T(k, to), 
effectively cuts the power on small scales and makes (p~3|) finite. In addition, today 
T(k,t ) becomes already negligible for scales k <C k e ; this is why the upper limit 
of integration in (|1"3D can be replaced by infinity. 

The ratio between Oij(to) an d ^_h"(^o) often used (see e.g. @, 



a H (t ) w 5 6 H {to), 
when simply substituted into (01), would yield the result 



a 2 H (tk) « 25 S 2 H (t ) 



~M H (t )~ 


n-1 
3 


~M H {t eq y 


_M H (t eq )_ 




_M H (t k )_ 



n-1 
2 



(14) 



(15) 



Eq. fll5|) is a reasonable first guess, but too inaccurate if one is willing to make 
precise predictions. According to (|9|) and (PH), one must instead use for t fc <C t eg 
the expression 



100 
~8T 

200 

9V 



« 2 (fc) S 2 H (t ) 



a 2 {k) kl $ 2 (k ,t ) 



' M H (t ) ' 


n-1 
3 


~M H (t eq y 


M H (t eq )_ 




_M H (t k )_ 



n-1 

2 



" Aftf (t ) " 


3 


"M^(t eg )" 


M H (t eq )_ 




_M H (t k ) _ 



n-1 
2 



(16) 



Eqs. (0), (0); an d (pTT|) are the main results of this section. Eq. (|15D corre- 
sponds, in our notation, to a{k) = a(k ) = 5 x 9/10 = 4.5. It is not correct 
as a(k) is both model and scale dependent, since the transfer function depends 
on the cosmological parameters and the time t k . Even for a given cosmological 
model for which (0) holds, it would certainly be incorrect at much earlier times 
tfc, according to (|Tl~D . This can result in a large discrepancy in the computation 
of (3(Mh)- It is recognized that the problem is much deeper than just a clever 
choice of the window function. Actually, the window function corresponding to 
the numerical results obtained in [T(| is the top-hat window function. If one uses 
8min ~ 0.7 as found by these authors, one must also use the top-hat window 
function in order to be consistent with their numerical results. 

It is tacitly assumed in (|16D, (as in (^j), that the PBHs form before the 
time of equality, t k t eq , and that the universe is first radiation-dominated 
and instantaneously becomes matter-dominated at t eq (with A = 0). A more 
complicated evolution of the scale factor a(t) will give a different expression. The 
observational input is the present density contrast <5#(to) on the Hubble-radius 
scale which is found using the CMB anisotropy data. 

What are the scales for which the discrepancy between ( jl6|) and (|15|) is most 
important? To get an idea, we consider the earliest formation times t k ~ t ke1 
right after inflation. This gives 



a 2 (k e ) 



x 1 



~ 2 W? 



TH 



'x) dx ss 0.21 



(17) 
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with 

0.202 < a 2 (k e ) < 0.218 for 1.3 > n > 1 . (18) 

Clearly there is a big difference for the earliest formation time t k ~ t^. A 
first rough estimate of a{k) corresponding to M H « 10 15 g gives (^) 2 a 2 (M H rs 
10 15 g) ^ 2, substantially smaller than 25. This estimate follows by taking into 
account that for modes outside the horizon one has \S^\ oc a 2 , while for modes 
within the horizon \5k\ is oscillating. . 

Interestingly, as seen from (||), the change implied by a smaller an{tk) at 
earlier times, amplifies the change induced by a larger 5 min - one is led to a 
smaller f3(M). Therefore, constraints on n max from PBH formation both on mass 
scales M < 10 15 g (these concern PBHs that have already evaporated by today) as 
on scales M > 10 15 g (see below) will be less constraining if one uses ( |TED instead 
of flip]). We note finally that an effective quantum-to-classical transition takes 
place already for PBHs with masses M -C 10 15 g , but not for PBHs with M ~ lg 
corresponding to k ~ k e ]5IJ (with the exception of PBHs possibly produced 
during preheating, see e.g. || [H|). At these formation times, the growing and 
decaying modes of the fluctuations that reenter the Hubble radius are still of the 
same order (no squeezing). Therefore, the classical approach must be amended. 
Anyway, a correct expression of the power spectrum must then take the decaying 



mode into account 21 



2.3 Inflationary spectra with a scale 

It is also our intention to investigate what happens when the assumption of a 
scale-free spectrum, tacitly or implicitly done in the literature, is dropped. We 
therefore consider spectra with a characteristic scale for which ([5]) and therefore 
(|16|) do not apply anymore. For each specific spectrum we have to compute Otf(tfc) 
numerically, but it is possible to give a rough estimate of what will happen if we 
consider some fiducial cases. Consider for example the case of a spectrum with a 
pure step in anitk) at k — k s , with tk a < t eq , where the ratio of anitk) on large 
scales to that on small scales is given by p. In this case, we have 



OHitk) = 77 oc(k) 8 H (t ) 

y 



M H (to) 
_M H (t. 



eg J J 



M H (t eq 
_M H (t k )_ 



X 



1 

p~ 



for k < k s 
for k > ke 



(19) 

We emphasize that, as explained in Sect. 2, this does not correspond to a 
pure step in the primordial power spectrum Sn{tk)- But with the approximation 
of a constant a for all PBH masses greater than 10 15 g, this is certainly the 
simplest extension to ([U|) that has a characteristic scale. Here we shall restrict 
ourselves to this case. In a future paper we shall then consider models with 
a characteristic scale which are motived by an underlying Lagrangian, or even 
features in the inflaton potential, for example a particular primordial fluctuations 
spectrum which derives from a jump in the inflaton potential derivative. The 
resulting fluctuations spectrum for such a potential is of a universal form and an 
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exact analytical expression has been derived by Starobinsky [Tj| . This expression 
depends (in addition to the overall normalization) on two parameters p and k s , 
which have a similar meaning to those in the toy model considered here. 

In our future work we will also take into account the effect of a cosmological 
constant corresponding to ~ 0.7 at present, since this seems to be strongly 
supported by recent observations and makes our current comprehension of the 
universe converge into a coherent picture. 



3 Observational constraints 

There are various limits on the initial mass fraction (3{M) of PBHs 0, most of 
them being related to the effects of Hawking radiation. Here, we restrict ourselves 
to the gravitational constraint 

QpBHfl = 2 &PBH,eq ^m,0 < 1 j (20) 

which states that the present PBH mass density must not exceed the present 
density of the universe. In (|2"0"D, Q m ,o refers to all dustlike matter today including 
a possible contribution from PBHs. In order to relate ( p0|) to the initial PBH 
abundance, one needs in principle to know the entire history of the universe right 
from the time when a mass scale M crosses into the horizon. Then one can 
compare the evolution of the PBH density ppbh (which goes like ppbh oc a~ 3 ) 
to the total density p to t- In a radiation dominated universe p to t oc a" 4 , so a small 
initial fractional PBH density can grow very large until today. We take again 
the evolution of the scale factor corresponding to fllEp . We follow steps similar 
to 0, but take into account all ^-factors, i.e., use the constancy of the entropy 
S = g*s{T)a 3 T 3 . This leads to the limit 

(3{M H ) < 1.57 x 1(T 17 (J£§^ 2 h 2 , (21) 

where we have used t eg = 4.36 x 10 10 s (Vt m flh 2 )~ 2 x (^) 6 f22|. Strictly speaking, 
the behaviour of the scale factor is slightly more sophisticated, but we leave that 
aside here. Anyway, as we will show in Fig.l, a small change in fl2~l|) does not 
matter. Finally, since black holes of initial mass Mppn < 5 x 10 14 g will have 
evaporated by the present day, the tightest constraint is obtained for M# ~ 10 15 

g- 

3.1 Scale-free spectrum 

In order to compare the constraints on the spectral index with earlier results we 
first consider a scale- free spectrum, n = const, on all scales. We further make the 
approximation a(k) w a(ko) ~ 4.5 in order to compare our results with previous 
work. As we have seen in Sect. II, this is not correct. The exact numbers can 
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n 

1.24 1.26 1.28 1.30 1.32 1.34 1.36 1.38 1.40 




lg/9 



Figure 1: The quantity (3 = /3(M = 10 15 g) is shown as a function of the spectral 
index n, for several values of a = a(M = 10 15 g). The straight lines represent the 
gravitational constraint and, for the purpose of illustration, the same con- 
straint weakened and strengthened by a factor of 100, respectively. The resulting 
maximally allowed value for n clearly depends only weakly on the precise value 
of the prefactor of the gravitational constraint (0), whereas it depends rather 
sensitively on the value of a. The value usually used in the literature, a ~ 4.5, 
overestimates the initial PBH abundance significantly and thus leads to much 
stronger constraints on n (i.e. n < 1.27) than would be expected for the more 
realistic choice of a ~ 1.5 (see section |2.2|) , which results in n < 1.32. 

only be calculated numerically. This is beyond the scope of our paper, but we 
shall get an idea about the changes by displaying the dependence of the obtained 
constraints on the used value for a(k), see Fig. [I]. 

The data obtained from COBE can be used |7j to normalize the spectrum (H) 

to 

°H{M Hfl ~ 10 56 g) = 9.5 x 10" 5 . (22) 

The above numbers are understood as to include the factor 10/9. In this way, 
the comparison with earlier results is straightforward. With the following setting 
of the parameters, 

h := 0.5, S min := -, (23) 

and by using equations ([5|) and (0), evaluated at M H = 10 15 g, one gets for the 
spectrum (^j the following upper bound on the spectral index: 

n < 1.27. (24) 

This result differs from n < 1.31 and is much closer to the constraints they 
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found due to the evaporation of PBHs, n < 1.24. This justifies our approach 
of only looking at the gravitational constraint - particularly since the details of 
black hole evaporation are still quite speculative. As mentioned above, the real 
constraints (i.e., taking into account the full expression ( |TT| ) for a) are expected 
to be somewhat weaker, cf. Fig. |]. 

We also had a look at the dependence of the result on the particular choice of 
the parameters (p3|) . What we found was a weak dependence on h and M H (i.e. 
the minimal mass to which the gravitational constraint fl2"T| ) applies): varying h 
from 0.1 to 1 as well as varying Mh from 10 14 g to 10 16 g changes the constraint 
on n by only about 0.01. 

The motivation for varying Mh is twofold: first, it amounts to the fact that 
the details of black hole evaporation are not yet fully understood, so the minimal 
mass of PBHs that have not evaporated by the present day is afflicted with some 
(small) uncertainty; second, if the mass of the PBH formed is not equal to the 
horizon mass Mh, but somewhat smaller, i. e. Mpbh — cMh, e < 1, the minimal 
horizon mass to which the gravitational constraint fl21~|) applies, is changed to 
Mh = t~ l ■ 10 15 g. Semianalytical considerations, for example, require e = 72 M. 
In the radiation dominated case (7 = |) this would weaken the constraint (|24"D 
by only about 0.005. But even e = 0.001 weakens it by only slightly more than 
0.02! 

On the other hand, we have found, as expected, the dependence on the value 
of S m i n to be rather strong: for 5 miri = 0.2 one gets n < 1.24 and for <5 min = 0.7 
one gets n < 1.30. The latter point is rather important, since 5 m i n = § (as 
assumed above) relies on semianalytical arguments, whereas numerical analysis 
seems to support 5 min ~ 0.7 [16|]. Fig. || shows the dependence of n max on 6 min . 

3.2 Scale-free spectrum with a superimposed step 

Next we consider the case of a simple extension to the scale-free spectrum, i.e. a 
power-law spectrum with a step as described by fliTf). Again, we use the approx- 
imation a(k) m a(ko) « 4.5 in order to compare our results with the previous 
subsection. We restrict ourselves to the case p < 1, i.e. more power on smaller 
scales. The opposite case leads to interpretational problems for the quantity 
(3{M): For p > 1 it is no longer a monotonically falling function of mass and 
therefore cannot be interpreted anymore as the fractional mass density of PBHs 
of mass larger than M. For such a case, the Press-Schechter formalism would 
have to be modified. But p < 1 is the physically interesting case for us, since it 
would produce more instead of less PBHs on smaller scales. 

In contrast to the scale-free case we now have three parameters - n, p and 
k s (or Mn(tk a ), respectively) - to be determined. One would therefore generally 
expect that the constraint ([H|) on (3{M) gives a plane in this parameter space. 
Since this constraint is a rising and (3{M) a falling function of mass, the strongest 
constraint will as before come from M ~ 10 15 g. But as long as M H (t ks ) > 10 15 
g, P(M ~ 10 15 g) does not depend on the choice of k s , so we get a functional 
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Figure 2: Dependence of the constraint on the spectral index, n max , on 5 min for 
the case of a scale-free spectrum, with h = 0.5, Mh ~ "'' ' 
n is clearly weakened for larger values of S min . 



10 g. The constraint on 



dependence between n max and p only. In order to illustrate this point, Fig. |3| 
shows (3{M) together with the gravitational constraint (21). 

Applying the same calculations as described in the previous paragraph, one 
gets n max (p) as depicted in Fig. f|. For p — 1 the result of the previous subsection 
is recovered (i.e. n max = 1.27), but for lower values of p the constraint on n 
may be considerably strengthened. For example, p = 0.01 requires n < 1.06 ! 
Of course this result can also be used the other way around: if n > 1.06 were 
found by some experimental measure, the gravitational constraint on the PBH 
abundance would require p > 0.01. 

The analysis sketched above does not tell us anything about the scale k s at 
which the step occurs. It would nevertheless be interesting to learn something 
about k s . This could arise, for example, from observations of MACHOs, hinting 
at the presence of a feature in (3{M) in the mass range M > 10 15 g. 



4 Conclusions 

We have considered two aspects of the formation of PBHs. On the one hand, we 
have derived a more precise expression for the PBH abundance. It corrects the 
usual formula used, and could in principle substantially improve the accuracy of 
the obtained results. The amount of correction implied would need numerical 
calculation and will be presented elsewhere. Our results concerning this part of 
the present work are summarized in (|16"D, (|T0"|), and (|TT|). On the other hand, 
we have considered simple models where the spectrum is not scale-free and have 
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compared it to the scale- free case. In the latter case we have found a somewhat 
stronger constraint than in [0] and have discussed the dependence of the result 
on the value of S m i n . In the former case, we have found that the presence of a 
step in the mass variance can drastically strengthen the constraints on n, giving 
in turn a tool to constrain the height of the step. This study can be extended 
to models with an underlying miscroscopic Lagrangian or phenomenologically 
motivated features in the potential. We plan to address these issues in a future 
publication. 
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Figure 3: This figure shows /3(M) for a step spectrum as in fll9|), with (arbitrarily 
chosen) n = 1.27, Mh^J = 10 17 g an d P = 0-5- The straight line is the 
gravitational constraint ( pi]) , which applies only for M > 10 15 g. The figure 
illustrates that the gravitational constraint is to be evaluated at 10 15 g and that 
the result does not depend on the choice of k s (as long as M H (tk s ) > 10 15 g). 
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Figure 4: Dependence of n max on p for the step spectrum (|E|). For p = 1 
the result for the scale-free case is recovered. For p < 1 the constraint on n is 
considerably strengthened. 
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